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The rate of quasi-stationary, two-dimensional magnetic reconnection is calculated in the frame- 
work of incompressible Hall magnetohydrodynamics (MHD), which includes the Hall and electron 
pressure terms in the Ohm's law. The Hall-MHD equations are solved in a local region across the 
reconnection electron layer, including only the upstream region and the layer center. In the case 
when the ion inertial length di is larger than the Sweet-Parker reconnection layer thickness, the 
dimensionless reconnection rate is found to be independent of the electrical resistivity and equal to 
di/L, where L is the scale length of the external magnetic field in the upstream region outside the 
electron layer, and the ion layer thickness is found to be di. 

PACS numbers: 52.35.Vd, 94.30.cp, 96.60.lv, 52.30.Ex 



INTRODUCTION 

Magnetic reconnection is a fundamental process of 
breaking and topological rearrangement of magnetic field 
lines in magnetized plasmas. Reconnection converts mag- 
netic energy into kinetic and thermal energy and is be- 
lieved to be responsible for many phenomena observed 
in the laboratory and cosmic plasmas [l|, 0j- Because 
electrical resistivity is very low in hot plasmas, mag- 
netic reconnection due to resistive dissipation of mag- 
netic field is typically a very slow process [EUS]. Re- 
connection can become much faster in the case when 
resistivity is anomalously high due to local plasma in- 
stabilities d, 0, [H • Another possibility is fast reconnec- 
tion made possible because of two- fluid plasma effects [l| , 
which require a two-fluid magnetohydrodynamics (MHD) 
description of plasma. In the limit of zero electron-to-ion 
mass ratio, two-fluid MHD equations simplify and reduce 
to Hall-MHD equations. The later include the Hall and 
electron pressure terms in the Ohm's law, in addition 
to the resistivity term present in single-fluid MHD. Be- 
cause of its relative simplicity, Hall-MHD description of 
plasma has been extensively used in numerical simula- 
tions of magnetic reconnection [(| @, H, 0- However, to 
the best of our knowledge, a simple self-consistent an- 
alytical model of Hall reconnection, similar to the clas- 
sical Sweet-Parker model of resistive reconnection, has 
not been constructed. In this paper we consider Hall- 
MHD equations and present a theoretical model of Hall 
reconnection. Full two-fluid calculations for magnetic re- 
connection will be considered elsewhere. The analytical 
derivations of this paper are similar to the derivations 
done by Malyshkin, Linde & Kulsrud || for the case of 
resistive single-fluid MHD reconnection. 



HALL-MHD EQUATIONS 

For simplicity and brevity, we use physical units in 
which the speed of light c and four times tt are replaced by 



unity, c = 1 and 4tt = 1 . To rewrite our equations in the 
CGS units, one needs to make the following substitutions: 
magnetic field B — > B/\/47r, electric field E — > cE/V^r, 
electric current j — > V^Lrj/c, electrical resistivity 7/ — > 
r]c 2 /4tt, and the proton electric charge e — > y/Ane/c. 

We assume the plasma is non-relativistic, with both 
phase and physical velocities much smaller than the speed 
of light. We neglect electron inertia for the Hall-MHD 
description of plasma. The generalized Ohm's law is 

E = - V x B + r/j + — j x B — • P e , (1) 

pe pe 

where is the ion mass, p is the plasma density, V is 
the plasma velocity (equal to the ion velocity), P e is the 
tensor of the electron pressure. The first two terms on the 
right-hand side of Eq. {1} are the single-fluid MHD terms, 
the third and fourth terms are the Hall and electron pres- 
sure terms. The equation of plasma motion is [Io| 

P^+p(V-V)V = -V-P+jxB, (2) 

where P is the tensor of the total pressure (equal to the 
sum of the electron and ion pressure tensors), and we 
neglect plasma viscosity. Equation ([2]) appears exactly 
the same as in the case of single- fluid MHD. Note that 
V • B = 0, and, for non-relativistic plasma, V • j = 0. 

We consider Hall magnetic reconnection in the classi- 
cal two-dimensional Sweet-Parker-Petschek reconnection 
layer, shown in Fig. [TJ The layer is in the x-y plane with 
the x- and y-axes being perpendicular to and along the 
layer respectively. All d/dz derivatives are zero. The 
thickness of the reconnection layer is 25, which is defined 
as the thickness of the out-of-plane current (j z ) profile 
across the layer. Note that 28 is approximately equal to 
the electron layer thickness, while the ion layer thickness 
2A can be much larger. Velocity Vl n is the plasma inflow 
velocity in the upstream region at point M, outside the 
electron layer. The magnetic field B m at point M is in the 
y-direction. The out-of-plane field B z is assumed to have 
a quadrupole structure (see Fig.[T]), in agreement with nu- 
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FIG. 1: Geometrical configuration of the reconnection layer. 
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(a nearly uniform "guide" field component of B z is taken 
to be zero). The reconnection layer is assumed to have 
a point symmetry with respect to its geometric center 
point O in Fig. [T] and reflection symmetries with respect 
to the axes x and y. Thus, the x-, y- and z-components 
of V, B and j have the following simple symmetries: 
V x (±x,Ty) = ±V x {x,y), V v (±x,Tv) = TV y (x,y), 
V z (±x,Ty) = V z (x,y), B x (±x,Ty) = TB x (x,y), 
B y (±x,Ty) = ±B y (x,y), B z (±x,Ty) = -B x {x,y), 
jx(±x,Ty) = ±j x (x,y), j y (±x,Ty) = Tj y (x,y) and 
j z (±x, =Fy) = j z (x, y). In the derivations presented below 
we will extensively use these symmetries. 



SOLUTION FOR HALL RECONNECTION 

Now let us list four assumptions that we make in this 
study. First, we assume that the plasma flow is incom- 
pressible inside the reconnection layer, p = constant. 
Second, we assume that the electrical resistivity rj is con- 
stant and very small, so that the Lundquist number is 
very large. Third, we assume that the reconnection pro- 
cess is slow and quasi-stationary, so that we can neglect 
all time-derivatives d/dt in Eqs. (HI)-© and below. This 
assumption is satisfied if the reconnection rate is slow, 
E z *C VA-Bm, and there are no plasma instabilities in 
the reconnection layer. Fourth, we assume that the elec- 
tron and ion pressure tensors are isotropic, so that the 
pressure tensors in Eqs. |T]) and ^ are scalars. 

Using Ampere's law and neglecting the displacement 
current, we find the x- and y-components of the current 
to be j x = d y B z and j y — ~d x B z . Here and below we 
use convenient notations d x = d/dx, d y = d/dy, d yy = 
d 2 /dy 2 , d^yy = d 3 /dxdy 2 , and etc. The z-component of 
the current at the reconnection layer central point O is 



jo = (jz)o = (d x B y - d y B x ) o « B m /5, 



(3) 



where we use estimates (d y B x ) -C {d x B v ) and 
(d x B y ) w B m /S at point O. 

Next, neglecting the time derivative in Eq. ©, the 
equation for acceleration of plasma along the reconnec- 
tion layer, in the y-direction, is p(V ■ V)Vy = —d y P + 



jzB x — j x B z . We calculate the first order partial deriva- 
tive d/dy of this equation at point O and obtain 



p(dyVy) 2 = -{dl y P) +J {dyB x ) 

~ —B m (d yy B v ) m + j (9 l 
= 2Bl/L 2 + 3o (d y B x ) . 



(4) 



Here, we used the fact that the pressure term is 
(dyyP)o ~ (dy y By/2) m - (d yy B z / 2) = B m (dy y B y ) m < 
0. Thus, the drop of pressure P along the layer is equal 
to the magnetic pressure drop of the parallel field compo- 
nent outside the layer. This result follows from the force 
balance condition for the plasma across the reconnection 
layer (in analogy with the Sweet-Parker derivations for 
slowly inflowing plasma), and its rigorous proof can be 
found in The last expression in Eq. (0} is obtained 
by defining the magnetic field external scale as 



L 



-2B m /(dy y B y ) m 



(5) 



This is the scale of magnetic field just outside the recon- 
nection electron layer (at point M) and can be interpreted 
as the length of the layer. The {d y V y ) derivative on the 
left-hand-side of Eq. ([4| can be estimated from plasma 
incompressibility condition at the O-point, 



{dyV y )o = ~{d x V x ) w Vi n /8, 



(6) 



where we use an estimate {d x V x ) ~ —Vi n /S. 

Next, the Faraday's law V x E = —dH/dt for the x- 
and y-components of a quasi-stationary magnetic field 
in two dimensions is dE z /dy = —dB x /dt «j and 
dE z /dx — dBy/dt w 0. Therefore, the electric field z- 
component E z is constant in space, and from generalized 
Ohm's law (TTJ) we obtain 

constant in space Rj E z = —V x B y + V y B x + rjj z 

+ (m p /pe)(j x By - j y B x ). (7) 

Now, we use Eq. to calculate E z at points O and M 
(see Fig. [IJ. At point O we have E z — r)j . At point M 
we have E z = V in B m + r/(j z ) m + (m p / pe)B m (j x ) m Rj 
V in B m + rj(j z ) m + {m p /pe)B m (d 2 y B z ) S 1 where we use 
an estimate (j x ) m Rj {d x j x ) S = (d 2 y B z ) 5. Equating 
these two expressions for E z at points O and M and ne- 
glecting the resistive term outside the reconnection layer 
at point M, rj(j z ) m Rj r/B m / L <C rjj , we obtain 



Vjo ~ V ln B m + (m p / pe)B m (d 2 B Z ) S. 



(8) 



Next, we calculate the second order partial derivative 
d 2 /dy 2 of Eq. (0) at point O. We have 

Rj 2(dyVy) (dyB x ) + V (d 2 yyJz ) 

-2(m p / pe)(d y jy) {dyB x ) 
Rj 2{d y V y ) (dyB x ) -r 1 {2j /L 2 ) 

+2(m p /pe)(d 2 xy B z ) (d y B x ) . (9) 
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Here, to obtain the final expression, we use the fact that 
the y-scale of the current j z , to a factor of order unity, 
is about the same as the y-scale of the outside magnetic 
field, j- x {d 2 yy j z ) w B-^{dl y B y ) m = -2/L 2 . This result 
can be understood by taking the d 2 /dy 2 partial deriva- 
tive of the Ampere's law equation j a sa B m /S, see Eq. ©, 
while keeping S constant because the partial derivative in 
y is to be taken at constant x = 6 (for details see Q). 

Next, we use the z-component of the Faraday's law. 
We have pa —dB z /dt = d x E y — d y E x . We calculate the 
d 2 jdxdy partial derivative of this equation at point O 
and use Ohm's law ([1]) for E x and E y . After tedious but 
straightforward derivations, we obtain 

pa -r) [{d xyxx B z ) + (d xyyy B z ) ] 

- [{d y B x )M x V z ) + (d x B y ) (d 2 y V z ) ] 

+ (m p /pe) [{d y B x ) (dl x j z ) + (d x B y ) (d 2 yyJz ) ] 

w 2 V (d 2 xy B z ) /5 2 

-2(m p /pe) [U /d 2 )(d y B x ) +j 2 /L 2 ] . (10) 

Here, to derive the final expression, we use (d xyxx B z ) pa 
-2(d 2 xy B z ) a /S 2 , (dt yyy B z ) « (dt yxx B z ) , (d 2 xxJz ) « 
-2j /S 2 , {d 2 yy j z ) w -2j /L 2 and (d x B y ) » j . We 
also use formula (d y B x ) (d 2 x V z ) + (d x B y ) (d 2 y V z ) = 
0. To prove it, let us consider the z-component of 
plasma motion equation l|2|). which is pV x (d x V z ) + 
ppy(9j,yi) = j x B y — j y B x . We calculate the d 2 /dx 2 
and d 2 /dy 2 derivatives of this equation at point O 
and obtain 2p(d x V x ) (d 2 x V z ) = 2(d x j x ) {d x B y ) and 
2p(d v V y ) (d 2 y V z ) = -2{d y j y ) {d y B x ) . Thus, we 
have {d 2 xx V z ) Q = -{d 2 xy B z ) {d x B y ) /p{d y V y ) < 0, 
(d 2 yy V z ) = (d 2 xy B z ) (d y B x ) /p(d y V y ) > and 
(d y B x ) (d 2 x V z ) + (d x B y ) (d 2 yy V z ) = 0. 

Now, we have six equations © , (gj) , © , ©, © 
and (|10|) . and we have six unknowns: j , <5, Vi n , {d y V y ) , 
(d y B x ) and (d xy B z ) . Thus, we can solve for all un- 
known quantities. For convenience of presentation, we 
express the solution in terms of the Alfven velocity Va — 
B m /yfp, the ion inertial length di = rriije^fp and the 
Lundquist number S = LVa/t]- The solution is 



jo f= 


a ^( 5 V3+2^/L 2 ) 1/2 , 


(11) 


6 p 


d L (5V3+ 2S 2 d 2 /L 2 ^ V2 , 


(12) 




a V3Va 2S 2 d 2 /L 2 ^ lf \ 


(13) 






(14) 


(d y B x ) p 


s ^i(5V3+2^/i 2 ) _1/2 , 


(15) 


(0 xy B z ^ o p 


d 2SdiB m / L 3 , 


(16) 






(17) 



The last equation gives the reconnection rate i£ z . In 
the limit 5* 3> 1 and di <ti L the reconnection rate 



is slow, _E 2 <C VAB m , and our assumption of a quasi- 
stationary reconnection process is self-consistent. Equa- 
tion ([H]) implies that the ions are accelerated up to 
approximately Alfven velocity Va along the reconnec- 
tion layer of length L. At the same time, the rate 
of electron acceleration along the layer at point O 

iS (d y V^) = (d y Vy) - {m i /pe)(dyiy) = (dyVy) + 

{V A di/B m )(d xv B z ) » (V A /L)(V3+2Sd 2 /L 2 ), where V e 
denotes the electron velocity. 

The ion layer thickness 2A can be estimated as follows. 
In the upstream region outside the ion layer at x = A 
ideal single-fluid MHD applies. Therefore, at x = A 
and y = the resistive and Hall terms in Eq. ([7]) can be 
neglected and E z pa VnB m , where Vr is the plasma inflow 
velocity outside the ion layer. Velocity Vr « E z /B m is 
called the reconnection velocity. It can also be estimated 
as Vr « {d y Vy)oA. Thus, 

We discuss the solution (jTTJ-fjTSJ) in the next section. 
Now, let us make an important remark. Our analyt- 
ical derivations involve an approximate solution of the 
Hall-MHD equations in the infinitesimal neighborhood 
of line OM across the reconnection electron layer (see 
Fig. [T]). All physical quantities, j a , 5, Vi n , (d y V y ) , 
{d y B, x ) , (d xy B z ) , B m and L, are defined either at 
point O (the layer center) or at point M (the upstream re- 
gion). In other words, all our derivations involve only the 
upstream region and the layer center, and we do not need 
to consider the downstream region for estimation of the 
reconnection rate. This "local" equations approach was 
first developed in [j| for single-fluid MHD reconnection 
with anomalous electrical resistivity, and this approach 
works for Hall-MHD reconnection as well. Note however, 
that we define the field scale L by Eq. ([5]), and its exact 
value, as well as the value of field B m in the upstream 
region, depend on the "global" solution of the Hall-MHD 
equations outside the reconnection layer. Both L and B m 
enter our model as parameters. Determination of their 
values requires numerical simulations of the global field 
configuration and is not considered here. 



DISCUSSION 

When di <C L/\/~S = Ssp (Ssp is the Sweet-Parker layer 
thickness), the solution (fnj) - (fT8]) reduces to the Sweet- 
Parker solution: j sa y/SB m /L, 5 pa A sa L/y/S = Ssp, 
V m « Va/VS, (dyVy)o w Va/L, (d y B x ) Pa B m /LVS, 
E z pa V A B m /VS and (d y V y e ) pa {d y V y ) pa V A /L. 

In the opposite limit, when d{ 3> L/\fS = <5sp and 
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reconnection is in a collisionless Hall regime, we have 



jo P 


s SdiB m /L 2 > %/S B ro /L, 


(19) 


S p 


a L 2 /5d,<L/\/5 = ( 5sp, 


(20) 


A p 


s 4» L/VS = 5 S p, 


(21) 


V5„ p 


s V A L/Sdt^V A /VS, 


(22) 


(8^)0 p 


a Va/L, 


(23) 


(dyB x ) P 


^ -B m / Sd{, 


(24) 


(0 x yB z ^ o P 


a 2SdiB m /L , 


(25) 


^ P 


a (di/L) Vk5 m » V A B m /^/S. 


(26) 



This run-away reconnection process can possibly operate 
in solar corona and Earth's magnetosphere. The depen- 
dence E z cx p~ x also supports a model of self-regulation 



The electron acceleration rate along the reconnec- 
tion layer at O-point is {d y Vy) ps Sd 2 V A /L 3 > 
(dyVy) . Fast electron outflow along the layer creates 
the quadrupole field B z [l8j]. From Eq. (|26|) we find 
that the rate of collisionless Hall reconnection, E z — 
(dj/L) V A B m , is independent of the electrical resistiv- 
ity 77 [12]. The reconnection velocity is Vr ~ E z /B m ps 
(di/L)V A . From Eq. (f2"Tj) we find A ps d t for the ion layer 
thickness, which is in agreement with experiment fl3l ]. 

In it noteworthy that in the absence of collisions, when 
77 — > 0, it is an anisotropic electron pressure that balances 
E z field at point O. In this case electrons are acceler- 
ated by E z field during time ~ L/V e T, while they are 
unmagnetized and are traveling with thermal speed V e T 
inside the electron layer of length L pil |20T|. As a result, 
the effective resistivity becomes 77 c ff ~ dlV e T/L, where 
d e = y/mim e / e-yfp is the electron inertial length (m e is 
the electron mass). Thus, pressure anisotropy and elec- 
tron inertia become important (they will be considered 
elsewhere) . 

In the end of this paper, let us briefly discuss possi- 
ble mechanisms of fast reconnection, which is indepen- 
dent of the macroscopic size of the reconnecting system. 
First, note that the collisionless Hall reconnection rate 
E z = (di/L)V A B m is high when the reconnection elec- 
tron layer length L is microscopically small and compa- 
rable to di. Although the value of L cannot be deter- 
mined in our model, numerical simulations find that L 
can indeed be much smaller than the macroscopic sys- 
tem size @, H Q EE E3]. For example, L can be 



as small as ps lOd, during quasi-stationary collisionless 
reconnection, resultin g in a very fast reconnection rate 
E z n0.1V A B m E~ 



Another possible mechanism of fast reconnection can 
result from the dependence of the reconnection rate on 
the density p and temperature T of the plasma. We 
have 77 cx T~ 3 / 2 [ljj, V A cx p' 1 / 2 , S cx p~ l / 2 T 3 / 2 , 
d 2 cx p^ 1 . Therefore, di/5sp cx (T/p) 3/i , and the col- 
lisionless Hall reconnection rate is E z cx p^ 1 . Thus, even 
if initially L> di, as reconnection proceeds, the plasma 
is heated up by Joule heating, the plasma temperature 
rises, the plasma density drops, and the reconnection be- 
comes more and more collisionless and faster and faster. 
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A note added after publication: After the publication 
of his paper, the author has became aware that Simakov 
and Chacon (2008) have independently obtained a Hall 
reconnection rate formula, which is the same as Eq. (17) 
up to numerical factors of order unity j2^|. Although 
the final result is basically the same, the two derivations 
are significantly different. The author has used a rigor- 
ous local analysis, in which the Hall-MHD equations are 
solved only in the upstream region and in the center of 
the reconnection layer. By contrast, Simakov and Cha- 
con (2008) have used a different approach that also con- 
siders the downstream region. In particular, they have 
assumed that the thickness of the ion layer is given by 
max{rfj, 5}, and obtained their reconnection rate formula 
as an approximate interpolation between the rates in the 
collisional and collisionless regimes. There are also minor 
differences in the problem formulation; that of Simakov 
and Chacon (2008) includes electron viscosity, which is 
neglected by the author, and neglects electron pressure, 
which is included by the author. Expression (di/ 'L)V A for 
the reconnection velocity in the collisionless regime was 
obtained earlier by Stanley W. H. Cowley (1985) 
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